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Thermal Finite Energy QCD sum rules for the vector current correlator are used to study quark- 
gluon deconfinement. Assuming p-meson saturation of the correlator in the hadronic sector, and 
the Operator Product Expansion in QCD, we obtain the temperature behavior of the resonance 
parameters (coupling, mass, and width), and of the leading vacuum condensates, as well as the 
perturbative QCD threshold in the complex squared energy plane. The results are consistent with 
quark-gluon deconfinement at a critical temperature Tc ~ 197 MeV. The temperature dependence 
of the p-meson width is of importance for current experiments on dimuon production in nuclear 
collisions. 
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I. INTRODUCTION 



The temperature behavior of the light-quark vector cur- 
rent correlator in the framework of thermal QCD sum 
rules (QCDSR) was first discussed in and later re- 
analyzed by others [2] mostly using Laplace transform 
sum rules. Inconsistent results from the use of these 
sum rules were first pointed out in A better 

understanding of the QCDSR method Q, both at zero 
and at finite temperature, and in particular the use of 
Cauchy's theorem in the complex energy plane to for- 
mulate quark-hadron duality, has led to a preference of 
Finite Energy QCDSR (FESR) over the Laplace trans- 
form counterparts. In fact, modern determinations of the 
QCD strong coupling @ and quark masses 01, as well 
as thermal properties of hadrons 's']-'^ are now mostly 
done in the framework of FESR. Regarding the latter, the 
emerging picture is as follows. A key parameter signaling 
quark-gluon deconfinement is the squared energy thresh- 
old, so{T), above which the hadronic spectral function 
is well approximated by perturbative QCD (PQCD), as 
first proposed in Explicit determinations of so{T) in 
several light- and heavy- light-quark systems 
find this parameter to decrease with increasing tempera- 
ture, vanishing at a critical value T = Tc interpreted as 
the deconfinement temperature. Two other hadronic pa- 
rameters, the coupling and the width, also behave as ex- 
pected from a deconfinement scenario, i.e. the coupling 
decreases and the width increases with increasing tem- 
perature. A monotonically increasing hadronic width, 
interpreted as resonance absorption in a hot plasma, was 



proposed long ago as a clear signal of deconfinement [l^l , 
and used to predict resonance broadening in dimuon pro- 
duction in high energy nuclear collisions . This effect 
was later confirmed by several experiments (l2j . An ex- 
ception to the above behaviour has been found for the 
heavy-heavy-quark states J/ip, Vc, and Xc, which appear 
to survive above Tc [l^ in agreement with lattice QCD 
results [3] . It should be recalled that the conceptual in- 
terpretation of a hadronic width at finite temperature is 
different from that at T = 0. While in the latter case the 
width is entirely related to decay into allowed hadronic 
channels, at T particles hadronically stable at T = 
develop a width due to the emergence of scattering chan- 
nels which modify the interaction rate. This has been 
demonstrated explicitly in hadronic models, e.g. the lin- 
ear sigma model [15| , chiral perturbation theory , and 
others, as well as in QCD sum rules [i] , where e.g. the 
pion and the nucleon develop a width which increases 
monotonically with temperature. In other words, an in- 
creasing width signals an increase in the interaction prob- 
ability, rather than only a decrease in a decay rate. 
Regarding the hadron mass, its temperature behavior is 
hardly of any interest in the context of the deconfine- 
ment transition. In fact, the mass is nothing but the real 
part of the pole of the hadron propagator in the complex 
squared energy plane. Whether the pole moves up or 
down with increasing T does not provide, in itself, any 
information regarding deconfinement. It is the imaginary 
part of the hadron propagator, i.e. the hadronic width, 
which signals deconfinement if it grows with increasing 
T. One could picture the extreme situation of the mass 
decreasing and vanishing at some T* . If this behavior is 
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not accompanied by an increasing width then the hadron 
would still be present in the spectral function at T = T* . 
Explicit determinations of the temperature dependence 
of hadron masses indicate a mild change with tempera- 
ture, increasing or decreasing slightly depending on the 
channel. 

In this paper we consider the first three thermal FESR 
for the vector current correlator to find the behavior of 
so{T) and the p-meson parameters, as well as the di- 
mension d = 4 gluon condensate, and the d = 6 four- 
quark condensate entering the Operator product Expan- 
sion (OPE) in QCD. This follows a recent analysis [17] of 
the axial- vector channel with a much improved hadronic 
sgectral function, involving not only the pion pole as in 
[l^, ^18^ but also the ai(1260) resonance. Since the pion 
pole in this channel is related to the quark condensate, a 
connection can be established between chiral-symmetry 
restoration and deconfinement. Given the absence of 
a pole in the vector channel, it is not a-priori obvious 
that results in this case should be similar to those in 
the axial-vector channel, where the thermal quark con- 
densate mostly drives the behaviour of so{T). Another 
difference between axial- vector and vector channels is the 
presence in the latter of a space-like cut in the complex 
energy plane, interpreted as due to the scattering of the 
vector current off pions in the hot plasma This term 
is of higher order (two-loop) in the axial-vector channel, 
and thus negligible On the other hand, 

PQCD is chiral symmetric (in the chiral limit), and so 
is the leading dimension d = 4 term in the OPE. Ad- 
ditional motivation for reexamining the thermal vector 
current correlator using FESR, and modern information, 
is to determine the T-dependence of the p-meson width 
entirely from QCD, as opposed to earlier work based on 
current algebra [111] . This is important in connection 
with current experiments on dimuon production in high 
energy nuclear collisions (T2j . 



where radiative corrections are currently known up to 
five-loop order, i.e. 0{ai). The QCD FESR rest on two 
pillars [3, (i) the operator product expansion (OPE) of 
current correlators at short distances beyond perturba- 
tion theory, and (ii) Cauchy's theorem in the complex 
squared energy plane, which relates the (hadronic) dis- 
continuity across the cut on the real semi-axis with the 
integral around a contour of radius \so\ where the OPE 
is expected to be valid. The latter is usually referred to 
as quark-hadron duality. This leads to the FESR 
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Imn(s) 



j^[l + 0{as)] (iV^l,2,...), 



HAD 



(3) 



where the leading order vacuum condensates in the chiral 
limit are the dimension d — -i gluon condensate 



^4(04) - -(a,G2), 

and the dimension d = 6 four-quark condensate 
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(4) 



^6(06 



((97mA^)') 



(5) 

The radiative corrections in Eq.(3) are known up to five- 
loop order, i.e. 0{a^), and they will be used at T = to 
normalize the FESR. 

Under the extreme approximation of vacuum saturation, 
the four-quark condensate can be related to the square of 
the quark condensate. However, there is no convincing 



II. QCD SUM RULES 



The light-quark vector current correlator at T = is 
defined as 




<0\T{V^{x)Vjm\0> 



= {-g^^.q^ + qt.qu)u{q^) , (i) 

where V^(a;) = i[: u{x)jfj.u{x) — d{x)^^d{x) :] is the 
(electric charge neutral) conserved vector current in the 
chiral limit, and — {^jJ,q) is the four-momentum car- 
ried by the current. The function n(q^) in PQCD is 
normalized as 

Imn(a = ^[l + 0(a,(g2))], (2) 
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FIG. 1: The normalized thermal behavior of the PQCD 
threshold in the vector channel (solid curve), for Tc = 
197 MeV, and in the axial- vector channel (dotted curve) from 
[l7t for the same Tc. 
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FIG. 2: The normalized thermal behavior of the p-meson lep- 
tonic decay constant, for Tc = 197 MeV. 



FIG. 3: The normalized thermal behaviour of the p-meson 
width, for Tc = 197 MeV. 



theoretical support for such an approximation. In fact, 
determinations of the vacuum condensates from data on 
hadronic decays of the r-lepton [l^, as well as e~e+ anni- 
hilation into hadrons f20l] , indicate strong deviations from 
vacuum saturation. Theoretical arguments from chiral 
perturbation theory also do not support this approxima- 
tion at next to next to leading order and at T = [2l| . 
An extension of this analysis to finite temperature indi- 
cates a breakdown of vacuum saturation except in the 
chiral limit . Since there are no gauge invariant oper- 
ators of dimension d = 2 in QCD, it is standard practice 
to assume they are not present in the OPE. This assump- 
tion is supported by results from data analyses (l9j-(20|. 
In the hadronic sector, assuming p-meson saturation of 
the spectral function, and a Brcit-Wigner resonance form 
gives 



1 , ,,11 
-Imll had(s) = --??■ 



M2r2 



(6) 



where [23] /p = 5 is the leptonic decay constant, Mp = 
0.776 GeV and 0.145 GeV are the p mass and width, 
respectively. This finite- width parameterization has been 
normalized such that the area under it equals the area 
under a zero width expression, i.e. lmIl\HADis) ~ 
JlMl5{s-Ml). 

A test of the FESR, Eq.©, with N = 1,2,3 can be 
performed by determining sq together with the gluon 
condensate and the four-quark condensate, and com- 
paring them with results from data analyses. Using 
the full PQCD information on n(s) to five-loop order, 
with a,(M2) = 0.338 ± 0.012 [6|, and Eq.® one finds 
So = 1.44 GeV^ 0^(6^) = 0.12 GeV and Ce{de) = 
—0.39 GeV^. These results are in reasonable agreement 
within errors with 19] and (20|. This is not surprising. 



as the condensate determinations based on experimen- 
tal data (T^-ji^] require similar values of so- The value 
y/so = 1.2 GeV validates the assumption of p-dominance, 
as the first radial excitation of the p is the p(1450) with 
a mass M ~ 1.5 GeV. 

The finite- width parametrization, Eq. (6), is clearly not 
unique. In order to test its impact on the results at 
T = we used instead of the standard p-propagator 
1/ [s- AP + iMT] the alternative 1/ [s- + i^/sV] . 
We find essentially the same solution to the FESR, at 
five-loop order, for sq and the d — A,6 vacuum con- 
densates provided the leptonic coupling fp is somewhat 
smaller, i.e. fp ~ 3. This will play no role at finite T, as 
we shall normalize all results to their T = values, i.e. 
we are essentially interested in finding the thermal be- 
havior of ratios. These ratios are hardly distinguishable 
from those using Eq.(6). 

The extension of the FESR program to finite temper- 
ature was first outlined in Field theory arguments in 
support of the validity of this extension were later given 
in 24]. At finite T there is an additional longitudinal 
structure in Eq.(l), but we shall consider the FESR for 
the transverse part. In the QCD sector one needs to 
restrict PQCD to the leading, one-loop level, as the ap- 
pearance of two scales in as{q^,T), i.e. Aqcd and Tc, 
remains an open problem (QCD sum rules approach Tc 
starting from T = 0, where PQCD is not valid). At this 
order there are two thermal contributions to the vector 
correlator, one in the time-like region {q^ > 0) and one in 
the space- like region (g^ < 0). In the static limit (q 0) 
these terms are 
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FIG. 4: The normalized thermal behaviour of the p-meson 
mass, for Tc = 197 MeV. 
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FIG. 5: The normalized thermal behaviour of the dimension 
d = 4 gluon condensate (solid curve), together with lattice 
QCD results [l|] (solid circles) for Tc = 197 MeV. 

for the time-like contribution, and 



Imn-(w,T) = 




= ^TH{u;'), (8) 

in the space-like region, where npiz) = 1/(1 -I- e^) is the 
Fermi thermal function, and the chiral limit was assumed. 
The vacuum condensates develop a T-dependence which 
can be obtained from the sum rules themselves, or by re- 
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FIG. 6: The normalized thermal behaviour of the dimension 
d = 6 four-quark condensate, for Tc = 197 MeV. 

sorting to lattice QCD (LQCD) determinations. A non- 
gauge invariant, non-zero dimension d — 2 term in the 
OPE only appears at high temperatures , beyond the 
domain being normally explored with QCDSR so that it 
can be safely neglected. 

In the hadronic sector the leptonic coupling, the mass, 
and the width of the p-meson entering Eq.(|ni) become 
temperature dependent. In addition, there is a hadronic 
contribution in the space-like region due to the coupling 
of the vector current to two pions in the thermal bath, 
and given by 

hmU-\HAD{uJ,T) = ^S{w^) ynB[^)dy, (9) 
where nsiz) = l/(e^ — 1) is the Bosc thermal function. 



III. RESULTS AND CONCLUSIONS 



Given that there are only three leading FESR, and 
a few more parameters, the following strategy has 
been adopted. First, results for the quark condensate 
from LQCD [2^ was used as a first approximation to 
Cg{Oq)(T) in vacuum saturation. Next, we took guid- 
ance from the temperature dependence of leptonic cou- 
plings, masses, and widths determined in other light- and 
heavy-light-quark hadronic channels. We then explored 
the six-dimensional parameter space using this informa- 
tion, seeking a solution to the three FESR that would 
make physical sense, i.e. that would agree with expec- 
tations as developed in other channels. This means a 
decreasing so{T), fp{T), Ci{Oi}{T), and Ce{06}{T), to- 
gether with an increasing T{T) with increasing T. A 
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scanning of the parameter space shows a welcome highly 
peaked structure with almost overlapping solutions in- 
volving two critical temperatures. One for deconfine- 
ment, Tc, pertaining to the vanishing of so(T), fp(T), 
and Ci{04)iT), and to the divergence of rp(T). And 
another temperature for the vanishing of the four-quark 
condensate, T*, which turns out to be essentially the 
same as that for the vanishing of the gluon condensate, 
and some 5% lower than Tc. This set of solutions re- 
quire the p-mass to decrease with increasing T, although 
it remains non-zero at T = Tc- However, since the res- 
onance width diverges at this temperature, the p-meson 
is no longer seen in the spectrum. The FESR cease to 
have solutions close to the deconfinement temperature, 
T/Tc ~ 0.90 — 0.95, as found in previous analyses of this 
channel [H-0, as well as in the axial- vector channel 
[T8| . The results are shown in Figs. 1-6, and correspond 
to the following analytical expressions 



rp(T) - 



rp(0) 



1 - {T/Tc)^ ' 
where a = 3, and = 197 MeV, 

Ce{d,){T) = cideKo) [i-iT/r;)"] 

with 6 = 8, and T* = 187 MeV, and 

Mp(T) = Mp{o) [1 - (T/r;,)^] , 



(10) 

(11) 

(12) 



where c — 10, and T^j — 222 MeV, constrained to satisfy 
> Tc. The slight 5% difference between Tc and T* 
is well within the accuracy of the method. A change 
in the values of the parameters a,b,c in Eqs. (ITUl) - P^ 
affects the behaviour of so{T), fp{T), and Ci{di){T). 
In order to retain the qualitative behaviour of the full 
six quantities, the parameters a, &, c are restricted to 
changes not greater than ±30%, ±50%, and ±30%, 
respectively. The temperature T* is rather tight, with a 
maximum allowed change of ± 3 MeV, while could 
vary in the range T;^ = 210 - 240 MeV. 



A fit to the results for the remaining three param- 
eters gives so(r)/so(0) = 1 - 0.5667 (T/Tj^-^s - 
4.347 (T/T,)68-4i, Ci{di){T)/Ci{di){Q) 
1 - 1.65 {T/Tcf-''^^ + 0.04967 {T/Tcf-™\ 
and /p(r)//p(0) = 1 - 0.3901 {T/TcY"-'' + 
0.04155 (T/rc)i-269, corresponding to Tc = 197 MeV. 
The behavior of so{T), Fig. 1, is somewhat similar 
to the recent result in the axial- vector channel [l^ . 
While the hadronic spectral function is very different in 
these two channels, PQCD is chiral-symmetric (in the 
chiral limit). This result is pointing to an approximate 
universality of the deconfinement transition in light- 
quark systems. The different behaviour close to the 
critical temperature can be traced to the contribution 
of the quark condensate (equivalently the pion decay 
constant) in the axial-vector channel, which is absent 
in the vector correlator (at d — 4 the term mq{qq) is 
negligible). The thermal width of the p- meson. Fig. 
3, exhibits a dramatic increase of roughly a factor 20 
near Tc. Its functional form, Eg. dTUl) . should be of use 
in current experiments measuring dimuon production 
in heavy ion collisions p^]. The result for the thermal 
gluon condensate. Fig. 5, is in good agreement with 
LQCD lj|, and the four-quark condensate, Fig. 6, is 
compatible with the behaviour of \{qq){T)\'^, albeit with 
a coefficient different from that in the vacuum saturation 
approximation. This coefficient, though, cancels out in 
the ratio. 
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